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ABSTRACT 



Context. Since solar-like oscillations were first detected in red-giant stars, the presence of non-radial oscillation modes has been de- 
bated. Spectroscopic line-profile analysis was used in the first attempt to perform mode identification, which revealed that non-radial 
modes are observable. Despite the fact that the presence of non-radial modes could be confirmed, the degree or azimuthal order could 
not be uniquely identified. Here we present an improvement to this first spectroscopic line-profile analysis. 

Aims. We aim to study line-profile variations of stochastically excited solar-like oscillations in four evolved stars to derive the az- 
imuthal order of the observed mode and the surface rotational frequency. 

Methods. Spectroscopic line-profile analysis is applied to cross-correlation functions, using the Fourier Parameter Fit method on the 
amplitude and phase distributions across the profiles. 

Results. For four evolved stars, /? Hydri (G2IV), 6 Ophiuchi (G9.51II), r] Serpentis (KOIII) and d Eridani (KOIV) the line-profile vari- 
ations reveal the azimuthal order of the oscillations with an accuracy of ± 1. Furthermore, our analysis reveals the projected rotational 
velocity and the inclination angle. From these parameters we obtain the surface rotational frequency. 

Conclusions. We conclude that line-profile variations of cross-correlation functions behave differently for different frequencies and 
that they provide additional information in terms of the surface rotational frequency and azimuthal order 
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1. Introduction 



stars: individual: e Ophiuchi, rj Serpentis, S Eridani, p Hydri - Techniques: spectroscopic - Line: 



After several previous claims, the first firm observational evi- 
de nce of solar-like oscilla tions in red-giant stars was presented 
by iFrandsen et al.l (l2002h for f Hydrae. This was followed 
by discoveries from two-site ground-based spectroscopic cam- 
paigns on e Ophiuchi (De Ridder et al. 2006) and rj Serpentis 
(IBarban et al. 200 4i). Based on theory of more luminous red 
giants dlDziembowski et al.l 1200 li) . the detected frequencies of 
these stars were interpreted as radial modes and the stars were 
modelled based on this assumption dHoudek & Goughl l2002t 
| De Ridder et al, 2006). 

iHekker et a D ll2006l) made the first time-series analysis of 
spectral line shape variations and attempted to perform spec- 
troscopic mode identification of the observed frequencies for 
three red-giant stars (e Ophiuchi, rj Serpentis and ^ Hydrae) and 
one subgiant (6 Eridani). They investigated line-profile varia- 
tions in the cross-correlation functions based on a pixel-by-pixel 
method, i.e., fitting a sinusoid at every velocity pixel across the 
profile. From the differences in the shape of the amplitude distri- 
bution across the profile for different frequencies they concluded 
that non-radial oscillations have to be present, although no defi- 
nite identification could be provided. 



Send offprint requests to: S. Hekker, 
email: saskia@bison.ph.bham.ac.uk 

* The software package FAMIAS developed in the framework of the 
FP6 European Coordination Action HELAS ( |http://www.helas-eu.org) 
has been used in this research. 



The CoRoT (Convection Rotation and planetary Transits) 
satellite performs photometry on parts of the sky for 150 con- 
secutive days with 32 / 512 second cadences in the centre and 
anti-centre direction of the galaxy. Among the observed stars 
are many red giants. With these observations the existence of 
non-radial oscillations in red-giant stars has been firmly proven 
(iDe Ridder et al.ll2009l) . This proof is based on the determina- 
tion of the harmonic degrees, using the asy mptotic relation for 
high-order low-degree modes (lTassoullll98 0). This relation pre- 
dicts that non-radial modes appear at regular intervals in the 
Fourier spectrum (large separation), with modes with harmonic 
degree i = I approximately half way in between the radial 
modes and modes with harmonic degree { - 2 close to the radial 
modes (small separation). Indeed, multiple 'ridges' are present 
in echelle diagrams (frequency vs. freq uency modulo la r ge sep - 
aration) of several stars presented by iDe Ridder et al.l (l2009l) . 
These results are confirmed by observations of red giants from 
the NASA Kepler sateUte, see .Bedding et all (120 Id) , and will be 
followed up. 



The asymptotic approximation is valid for p modes and can 
only be applied for giants and subgiants when the oscillation 
modes are trapped in the outer parts of the star. The trapping 
depends on the internal stellar structure. For stars in certain evo- 
lutionary states, it is inefficient, such that g modes and mixed 
modes are also observable, see e.g. Dupretetal. (2009). This 
degrades the regular pattern, and the asymptotic relation is less 
useful as a guide to perform identification of the mode degree {. 
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Table 1. Stellar parameters of e Ophiuchi, 77 Serpentis, 6 Eridani and j8 Hydri: Effective temperature {T^g) in Kelvin, 
radius (R) in R0, mass (M) in Mq, metallicity ([Fe/H]), surface gravity (log g) in c.g.s units, rotational velocity (i/sim), 
microturbulence (^micro) and macroturbulence (^macio) in kms"', parallax (n) in mas and the apparent magnitude in the V 
band (nii,). 



parameter 


e Ophiuchi 


;/ Serpentis 


6 Eridani 


P Hydri 


Teff [K] 


4970 ± 80^ 


4955 ± 80° 


5050 ± 100" 


5872 ± 44^ 


R[Rq] 


10.4 ± 0.45 


5.914 ± 0.094*+'' or 5.794 ± 0.097^ 


2.33 ± 0.03' 


1.814 ±0.017^ 


M[Mq] 


2.0" 


2.42±0.2^orl.72'"j5[|« 


1.1215' 


1.07 ±0.03^ 


[Fe/H] 


-0.07 ±0.11° 


-0.15 ±0.11° 


0.13 ±0.03'' 


-0.2* 


log g (c.g.s.) 


2.9 ±0.15° 


3.2 ±0.15° 


3.77 ±0.16'' 


3.952 ± 0.005^' 


us'mi [kms"'] 


3.5 ±0.5° 


0.44 ± 0.44° 


2.3 ±0.5" 


3.3 ±0.3" 


^micio [kms"'] 


1.5 ±0.1° 


1.3 ±0.1° 






^miicio [kms"'] 


3.55 ± 0.45° 


3.52 ± 0.45° 






n [mas] 


30.34 ± 0.79'' 


52.81 ±0.75'' 


110.58 ±0.88'' 


133.78 ±0.5'' 


m,. [mag] 


3.24 ± 0.02'' 


3.23 ± 0.02'' 


3.52 ± 0.02'' 


2.82 ± 0.02'' 



° [Hekker & Melendez (2007), ' Richic hi et aP j200^. " iKallinger et alJ ( l2008h . '' iPerrvman & ESAl fT997h. " iReiners & SchmittI 
( l2003l'). ^ Valenti & Fischer (2005), « Val enti & Fischeii j2005h isochrone stellar mass. '' ISantos et all ( 120041) ' iThevenin et alj ( |200C 
^ iNorth et al. ( 2007) , * Dravins et al. ( 199^ 



For non-radial solar-like oscillations it is expected that all 
azimuthal orders are excited. However, it is usually not possi- 
ble to distinguish between oscillation frequencies of modes with 
the same radial order (n), harmonic degree ({), and different az- 
imuthal order (w), due to the generally slow rotational velocity 
of evolved solar-like oscillators, the limited time span of ob- 
servations, and stochastic nature of the oscillations. Slow ro- 
tation induces only a very small frequency splitting, while the 
stochastic nature causes the oscillation frequencies to appear as 
Lorentzian profiles in the Fourier spectrum with a width depend- 
ing on the lifetime of the mode, i.e., shorter mode lifetimes im- 
ply broader profiles. We need long time series of data to obtain 
sub-juHz frequency resolution and to resolve the modes, which 
can have life times of the order of a few to hundreds of days 
(iDupret et al.l i2009). So, the frequency of a non-radial oscilla- 
tion mode is extracted from a broad profile consisting of several 
modes with different m values. Which of these modes is domi- 
nant depends on the inclination angle. In stars seen at high in- 
clinatio n angle ot ^ modes h ave higher visibility than m = 
modes tGizon & Solankil|2003h . Furthermore, for individual re- 
alisations, as we consider a single time series of observations, 
the stochastic interference with noise might boost or diminish a 
mode. Therefore, we are in principle able to observe any of the 
present azimuthal orders or combinations thereof. The latter will 
also behave different than m = modes due to the influence of 
m i^Q modes. 

Line-profile variations are mainly sensitive to the azimuthal 
order of a m ode and much less sensitive to its harmonic de- 
gree (see e.g. lAerts & E ver 2000). The determination of a non- 
zero azimuthal order would indicate that the oscillation mode is 
non-radial. This can be important in cases where the asymptotic 
approximation fails or provides ambiguous results, i.e., where 
different ridges in an echelle diagram can not unambiguously 
be identified. In case f is known from the asymptotic approx- 
imation the m value can put additional constraints on the in- 
ternal structure. Also, we can obtain the inclination angle for 
non-radial modes and the projected rotational velocity from line- 
profile analysis. These parameters allow to determine the stellar 
rotational frequency for non-radial oscillators. 

With this potential of the analysis of line-profile variations in 
mind, we have improve d the spectros copic line-profile analysis 
originally presented bv iHekker et alj (l2006l) (Sectio n 3). These 
improvements were necessary as Hekker et al.l (12006 ) were only 



able to match observations and theory qualitatively. In this work 
we perform a quantitative analysis. We use FAMIAS (Freq uency 
Analysis and Mode Identification for Asteroseismology, IZimal 
2008), a package of state-of-the art tools for the analysis of 
photometric and spectroscopic time-series data, to perform line- 
profile analysis (Section 4) on three stars, /3 Hydri (G2IV) 6 



ICarrieij 



Eridani (KOIV) and e Ophiuchi (G9.5III), for which previous 
mode identification was available (Bedding et ai] |200"7[ f 
I2OO2I : [Carrier et alJl2003l : iKalhnger et al. 2008). 

In Section 5, we analyse the observed line-profile variations 
for oscillation modes of /3 Hydri, 6 Eridani and e Ophuichi. For 
non-radial modes in these stars, as identified from the asymptotic 
approximation, we obtain the azimuthal orders and a constraint 
on the inclination angle. The latter combined with the projected 
rotational velocity gives an estimate of the stellar rotational fre- 
quency. 

The results of the line-profile analysis for /3 Hydri, 6 Eridani 
and e Ophiuchi gave us confidence that we can perform a quanti- 
tative comparison between line-profile variations from observa- 
tions and synthetic spectra, which provides the azimuthal order 
and for non-radial oscillators, the inclination angle and with that 
the surface rotational frequency. Therefore, we perform the same 
analysis for rj Serpentis (KOIII), for which no previous mode 
identification was available (Section 5). 



2. Observations 

For e Ophiuchi and 77 Serpentis we have spectra from the fiber- 
fed echelle spectrograph CORALIE mounted on the Swiss 1.2 m 
Euler telescope at La Silla (ESO, Chile) at our disposal. These 
were obtained during a two-site campaign using CORALIE and 
ELODIE, the fiber-fed spectrograph mounted on the French 
1.93 m telescope at Observatoire de Haute Provence, France, 
during the summer of 2003. The ELODIE spectra are available, 
but their cross-correlation profiles have a resolution of 1 kms ' 
which is coarse compared to the 0.1 kms"' resolution of the 
CORALIE spectra, and are hence not used for the present anal- 
ysis. The observations of 5 Eridani were taken with CORALIE 
during a twelve day campaign in November 2001. For /3 Hydri 
we use spectra obtained with the fiber-fed echelle spectrograph 
HARPS mounted on the 3.6 m telescope at La Silla (ESO, 
Chili), which were obtained during a two-site campaign using 
HARPS and UCLES, the echelle spectrograph mounted on the 
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3.9 m Anglo-Australian Telescope, Siding Spring Observatory, 
Australia, in September 2005. Spectra from UCLES are con- 
taminated with iodine absorption lines, used to obtain accurate 
radial velocity variations, and are therefore not useful for the 
present analysis. For a detailed description of the data we refer to 
the publications de scribing the observations of solar-like oscilla - 
tions in these stars, De Ridder et al. (2006): Barban et al.'('2004: 
[Carrier (2002): Carrier etal. (2003): Bedding et al. (2007), re- 
spectively. Stellar parameters of the four stars are listed in 
Table □ 

We note here that w e have data for f Hyd rae at our disposal, 
but as already shown bv lHekker et al.l (l2006l) the signal-to-noise 
ratio of the oscillations is too low for the present analysis and 
therefore we did not include ^ Hydrae in the current investiga- 
tion. 

3. Improved analysis 

iHekker et alj ( |2006|) used the INTER-TACOS (INTERpreter for 
the Treatment, the Analysis and the Correlation of Spectra) soft- 
ware package developed at Geneva Observatory (Baranne et al. 
11996 ) to compute cross-correlation functions (CCFs) of the ob- 
served spectra. The CCFs are computed using a mask contain- 
ing absorption lines of a KO star for e Ophiuchi, rj Serpentis and 
6 Eridani and have a resolution of 0.1 kms For the HARPS 
data of /3 Hy dri we use th e CCFs computed with the HARPS 
pipeline ( Ru pprecht et al.l 12004). which uses a mask with ab- 
sorption lines of a G2 star. These CCFs have a resolution of 0.5 
kms"'. 

Here we use these CCFs as a starting point of a line-profile 
analysis, but we apply three corrections before attempting to 
analyse the line-profile variations. These corrections appeared to 
be necessary, as effects not intrinsic to the star (instrumental ef- 
fects / changing weather conditions) cause additional changes in 
the line profiles on a night b y night basis. These e ffects hampered 
a quantitative analysis by .Hekker et al.l (l2006l) . The following 
corrections are applied to isolate the line-profile variations due 
to oscillations from the extrinsic ones: 

- The first correction concerns a continuum normalisation. 
Continuum on both sides of each CCF is selected and a 
linear polynomial is fitted through these continuum points. 
The full cross-correlation profile is then divided by this 
linear fit. 

- The second correction involves the isolation of the line- 
profile variations and removal of outliers. For each night of 
observations we calculate for each velocity bin the average 
CCF value and its standard deviation. In each bin the average 
CCF is subtracted from the observed continuum corrected 
CCF to isolate the variations. Outliers are identified as 
values larger than three times the standard deviation. CCFs 
with outliers in one or more velocity bins are fully discarded. 

- Finally, the flux variations in each velocity bin are imposed 
on the mean CCF of all observations. This is not strictly nec- 
essary to perform line-profile analysis, but we aim to use mo- 
ments to obtain the oscillation frequencies and therefore we 
need to have a line profile rather than the residual profile. 

The resulting average CCFs and isolated line-profile varia- 
tions for all four stars are shown in Fig. [T]together with grayscale 
plots of the isolated line-profile variations. 

Extensive tests have been performed to investigate the influ- 
ence of these corrections and removal of outliers, in particular to 



confirm that these indeed remove effects not intrinsic to the star 
and that these do not influence the intrinsic behaviour. These 
tests indicate that the results presented here are robust to degrad- 
ing the signal and to changes in the average profiles over night 
due to for instance changed weather conditions or instrumental 

effects. 

With FAMIAS (Zima 2008), we compute frequencies from 
the first moment (Aerts et al . 1992) for the new CCFs, similar 
to what was done by Hek ker et al.l (j2.006) for the uncorrected 
CCFs. In general we recover the same frequencies (or 1 day 
al iases) as Hekker et al. (2006) for e Op hiuchi, 77 Serpenti s and 
as' Bedding et al.l(l2007D : lCarrieil(l2002l) : ICarrier et al.l(l2003D for 
y6 Hydri and 6 Eridani, respectively. 

4. Line-profile analysis 

The line-profile an alysis is pe rformed with the Fourier Parameter 
Fit (EPF) method ( Zimal2006 ). The EPF method relies on a fit of 
the observational Fourier parameters across the line profile. For 
every detected pulsation frequency, the zero point, amplitude, 
and phase are computed for every velocity bin across the profile 
and compared with the same quantities from synthetic profiles. 
This procedure only relies on the velocity eigenfunction at the 
stel lar surface and is t herefore model independent (see Chapter 
6 of Aerts et al. (2009)). We use the standard sign convention for 
the azimuthal order m in FAMIAS, i.e., m > implies pro-grade 
modes. 

One has to keep in mind that this method is mainly sensi- 
tive to the azimuthal order (m) of the mode and much less to 
its harmonic degree (f), as modes with the same order but sim- 
ilar degree w ill give esse ntially the same amplitude and phase 
distributions (IZimall2006l) . 

For the present analysis we use the frequencies of the first 
moment of the new CCFs and stellar parameters as listed in 
Table [T] as input parameters. We first computed pulsationally 
independent parameters, i.e., equivalent width (EW), macrotur- 
bulent broadening (^macro), projected rotational velocity (i/sin;) 
and radial velocity (RV) and fix these. The inclination angle is a 
free parameter in all calculations. Then we use the EPF method 
to fit the amplitude and phase distributions either for modes for 
which we know the degree from previous determinations or for 
a grid of modes with degree Q < { < 4, and order -{ < m < £. 
Modes up to degree 4 were investigated as in spectroscopy the 
partial cancellation effect would theoretically allow this. For ex- 
amples of amplitude and phase distributions, see e.g., Fig.|2l 

Although we investigate stochastic oscillations, we did not 
include damping and excitation in the synthetic profiles, but take 
the stochastic effects into account in the interpretation of the re- 
sults. The main reason for this is that we know that damping and 
excitation c an cause asymmetri es in the amplitude and phase 
distribution (' Hekker et al.l 12006). but the exact nature o f these 
asymmetries varies per realisation. In their Fig. 10, Hekke r et al.l 
(2006) show examples of amplitude and phase distributions of 
simulated noise free line profiles with a two day damping time. 
For each of the modes £ = 0, 1,2 and m - - £ they show 10 
realisations. From these simulations it is clear that for different 
realisations of the stochastically damped and re-excited oscil- 
lations, the height of the amplitude distributions varies for all 
modes. In addition to that, the centroids of the distributions for 
modes with m + can be shifted from the laboratory wavelength. 
The asymmetries in the amplitude profiles are most pronounced 
for sectorial modes with £ - |m|, but are also present for other m 
+ modes. Moreover, from simulations of spectra with signal- 
to-noise ratios of 50, 100 and 150, we concluded that the asym- 
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Fig. 1. Average cross correlation functions (top), isolated line-profile variations after subtraction of the average CCF (middle, only 
one night is displayed for clarity) and grayscale plots of the isolated line-profile variations as a function of phase of the respective 
dominant frequencies, where darker regions indicate negative values, while lighter regions indicate positive values (bottom). From 
left to right fi Hydri, e Ophiuchi, rj Serpentis and 5 Eridani. 



metrics in the amplitude distribution are also present for modes 
with m - 0. In the interpretation of our results we take the asym- 
metries and shifts into account by shifting the computed ampli- 
tude distributions such that their central minima overlap with the 
central minimum of the observed distribution. The same shift is 
applied to the computed phase distributions. 

To identify the best fit to the data, values have been com- 
puted, as foreseen in FAMIAS. These are used to select the fits 
with below 1 . However, it turned out that a reasonably large 
number of profiles satisfies this requirement and that among 
these best fits the optimum m value can differ when investigat- 
ing only the central part of the line, or when the full line profile 
is taken into account. Therefore, we discuss the results of our 
line-profile modelling from visual inspection of the amplitude 
and phase distributions across the line pr ofiles as it has been 
don e previously in the l i teratii re, see e.g. iTelting et"al] (Il997h 
and 'Briquet et all (I2005L l2009t) for p Cephei, 6 Ophiuchi and 
HD 180642, respectively. 



5. Results 

The analysis of line-profile variations for solar-like oscillations 
in evolved stars as discussed above is applied first on stars 
for which mode identification has already been performed, i.e., 
jS Hydri, S Eridani and e Ophiuchi. 



Table 2. Global line profile parameters as computed with 
FAMIAS and used for the line profile analysis. Terminology 
similar to the one used in Table[T] 



parameter 


6 Ophiuchi 


rj Serpentis 


5 Eridani 


/3 Hydri 


EW [km s-'] 


3.0 


2.7 


2.8 


3.2 


V sin ( [km s"'] 


5.7 


5.6 


4.9 


4.3 


^macro [km S ] 


3.7 


3.4 


3.6 


2.6 



40 F 



30 r 

to 
o 




Inclination angle [°] 

Fig. 4. Histogram of the inclination angles of all computed mod- 
els for y6 Hydri in which is used as weight. See text for more 
details. 
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velocity [km/s] 



velocity [km/s] 



velocity [km/s] 





jllHt t^t^: 


iti|| 






Fig. 2. Amplitude (top) and phase (bottom) distributions as a function of velocity across the line profile of /3 Hydri for the three radial 
modes with frequencies obtained from the first moment (v): V(v> = 86.8 cd"' (1004.4 juHz) (left), V(v> = 73.0 cd"' (845.0 /iHz) 
(centre) and v^v) = 96.7 cd"' (1118.9 //Hz) (right). The mean radial velocity of the star is approximately 23.2 kms"'. The data is 
indicated with black dots and the (t, m) - (0,0) fit is indicated in red solid line (colours only in the online version). In cases where 
the observed distributions are not centred around the midpoint of the CCF, the fits are shifted (see Section 4 for further explanation). 






velocity [km/s] 



velocity [km/s] 



Fig.3. Same as Fig.|2]for two frequencies obtained foryS Hydri from (v): V(v> = 89.3 cd"' (1033.1 fjHz) (left), v<v> - 103.0 cd" 
(1192.3 juHz) (centre), for which the harmonic degrees have been determined previously to be ^' = 1. In these panels the (i,m) ■ 
(1,0) and (£,m) - (1, +1) fits are indicated with red solid and green dashed lines, respectively. 



5.1. /i Hydri 

For y6 Hydri the degrees of the modes with high est power have 
be en identified u s ing th e asymptotic relation dTassoull Il980t) 
by iBedding et aP (l2007h . From the HARPS data that we have 
at our disposal we obtain five frequencies in common with 
[Bedding et al. (2007), or one day aliases, for which we were able 
to an alyse the line p r ofile v ariations. For three common frequen- 
cies, IBedding et all (l2007h determined { - and for the other 
two modes £ =1 was found. We use these five frequencies for the 
line profile analysis. The results of the best fits to the observed 
line-profile variations obtained with the FPF method are shown 
in Figs. |2] and [3] We list the equivalent width, projected rota- 
tional velocity (v sin /) and macro turbulence ^macm for all four 



stars in Table |2] The v sin / we obtain here for the line profile of 
jS Hydri is shghtly higher than the independent determination by 
iReiners & Schmittl(l200l . 

Inspection of the amplitude and phase distributions reveals 
clearly that for the radial modes the fits to the observed line- 
profile variations are consistent within their errors (see Fig. |2]l. 
For the two modes with f = 1 we fitted -£ < m < { and plot 
the fits for the zonal mode and modes with m - \ and m - -I 
for oscillations with frequencies 1033.1 juHz and 1192.3 //Hz 
respectively (see left and centre panels ofFig.[3]l. Forv= 1033.1 
//Hz the fit with m = seems to be the best fit, while for v = 
1 192.3 //Hz the fit with m = -1 seems favourable over the m = 
fit. Nevertheless, for both frequencies the zonal and sectoral 
mode are consistent with the observations within the errors. 
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Fig. 6. Same as Fig.|4]for 6 Eridani. 



Because of the non-radial nature of the latter two modes, the 
line profile analysis depends on the inclination angle of the star. 
This inclination angle together with the projected rotational ve- 
locity leads to a surface rotational frequency (O). For the inclina- 
tion angle we take the interval defined by the weighted mean and 
standard deviation. The weight is defined as Wk — Xp/xl' where 
is th e of the best solution, see for more details Desmet et aH 
(l2009l) . For this star, we find an inclination of 55 + 17°. An his- 
togram of the weighted inclination is shown in Fig. |4] So for /3 
Hydri we find a surface rotational frequency ranging between 
3.6 yuHz and 5.5 yuHz for inclination angles between 72° and 38° 
and a projected rotational velocity of 4.3 km s"' as obtained with 
FAMIAS. 

5.2. 6 Eridani 

ICarried (|2002|) : [Carrier et al] (l2003h analysed the radial veloc- 
ity variations of 6 Eridani from CORALIE observations and 
found evidence for 15 oscillation modes whose harmonic de- 
grees could be identified using an echelle diagram. We obtain 
three frequencies in common (or 1 day aliases) in the variations 
of the first moment of the same data set, with enough signal-to- 
noise ratio to analyse the line-profile variations. The harmonic 
degrees of all three oscillation frequencies are known and used 
as input for the present line-profile analysis. The results of this 
analysis are listed in Table|2]and shown in Fig.|5] 

The most dominant mode we obtain from the first moment 
has a f requency o f 675.8 fuHz and was identified by I Carriej 
(l2002h : ICarrierer al. (2003) as a radial mode. Our fit of a radial 
mode, is indeed consistent with the data within errors. The fact 
that the 'centre' part of the amplitude distribution does not drop 
to is an indication of the presence of a non-radial mode. Indeed 
a fit with {£, m) = (1,-1) does provide a better fit to the amplitude 
distribution and the centre part of the phase distribution (see left 

panels of Fig. |5]l. 

The second frequency is identified by ICarrie 1 (I2002h : 

[Car rier et al. (2003) as an = 2 mode and from our present anal- 
ysis we see that (i, m) = (2,1) seems the most likely mode identi- 
fication (see centre panels of Fig.lH ). The third fre q uency is con- 
sistent with ^ = 1 as determined bv lCarrie^ (|2002|) : [Carrier et al.l 
(12003), but we are not able to discriminate between m = or m 
= -1 due to the lower signal to noise ratio, which mostly influ- 
ences the phase distribution of this frequency (see right panels 
ofFig.Ell. 



40 r 
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Fig. 8. Same as Fig.|4]for e Ophiuchi. 

Also for this star we investigated the weighted mean and 
standard deviation of all inclination angles of the computed syn- 
thetic profiles (see Fig.|6]for an histogram in which the weights 
are taken into account) to compute the surface rotational fre- 
quency. The inclination angle is 57 + 19°, which results in a 
range of surface rotational frequencies between 3.1 and 4.9 //Hz 
at inclination angles ranging from 76 to 38° and a projected ro- 
tational velocity of 4.9 km s"' . 

5.3. e Opiiiuclii 

iDe Ridder eTaTI (l2006l) and lBarbanet all (l2007h interpreted the 
frequencies observed in e Ophiuchi, obtained from radial veloc- 
ities derived from CORALIE and ELODIE and from photomet- 
ric space data (MOST) respectively, as rad ial modes. Based on 
the combined data set Kallin ger et al.l ('2008) interpreted the fre- 
quencies of e Ophiuchi as both radial and non-radial oscillation 
modes. They performed mode identification by comparing the 
observed frequencies with those of theoretical models to derive 
the degree of the modes. Hence this mode identification is model 
dependent. From the (same) CORALIE data we have at our dis- 
posal we obtain three frequencies of which only for one fre- 
quency the ha rmonic degree of the oscillation has already been 
identified by Kallinger et al.l (l2008l) . The results of the best fits 
we obtain with the EPF method are listed in Tableland shown 
in Fig. Q 

A first inspection of the amplitude and phase distribution im- 
mediately shows that m = -2 provides for all frequencies the 
least likely mode identification. 

For the 1 " frequency we find that both fits for m = as well 
as m = -1 are within the errors consistent with the amplitude 
distribution for the full CCF, where the m = -1 fit is superior in 
the centre of the line (see left top panel of Fig.|7]i. This is also the 
case in the centre of the phase distribution, where the m - -I fit 
matches the observations closely, while the m = fit lays outside 
the error bars. In the wings of the CCF the opposite is the case, 
i.e., the m = fit matches the observations and the m - -I fit 
falls outside the error bars. Thus from the line centre we would 
conclude that m- -1 and from the fine w ings that m = 0. Both 
are consistent with 'Kallinger et al.' ('2008') and only the l atter i s 
consistent with De Ridder et al. (2006) and Bar ban et al.l ( l2007h . 

The observed amplitude distribution of the 2'^ frequency 
shows a clear dip in the centre which is matched very closely 
by the m - Q fit, while the m = -1 fit lays just within the error 
in this region (centre top panel of Fig. |7]). The m = fit also 
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Fig. 5. Same as Fig. |2] for the significant frequencies of 6 Eridani obtained from (v): v<v> = 58.4 cd"' (675.8 fiHz) (left), V(v> = 
61.1 cd"' (706.9 yuHz) (centre) and V(v> = 57.5 cd"' (665.9 yuHz) (right). The mean radial velocity of the star is approximately 
-6.3 km s ' . The fits with m = 0, +1, + 2 are indicated with red solid, green dashed and blue dotted lines, respectively. 




Fig. 7. Same as Fig.|5]for the significant frequencies of e Ophiuchi obtained from (v); v<v> - 5.03 cd ' (58.2 fjHz) (left), V(v> = 
4.65 cd ' (53.8 //Hz) (centre) and V(v> = 6.46 cd"' (74.8 //Hz) (right). The mean radial velocity of the star is approximately 
-9.4 km s-'. 



matches the phase distribution best in both the centre and the 
wings, which implies this is the most likely mode identification. 

The asymmetry in the observed amplitude distribution of the 
3'''' frequency hampers a good fit, but in the line centre both m 
= and m - -I fall well within the error bars (right top panel 
of Fig. I?]). The centre part of the phase distribution is best fitted 
with the m = -1 fit while the line wings are again best fitted with 
an OT = fit. Thus here, the centre of the line favours a m = -1 
mode identification, but m -Qis also possible. 

For this star the preferred inclination angles, computed from 
the weighted mean and standard deviation, all lay in the range 
63 + 18°, which, together with the projected rotational velocity 
of 5.7 kms ', results in an surface rotational frequency ranging 
between 0.8 and 1.1 //Hz. 



5.4. 1] Serpentis 

After confirmation that we can compare the line-profile varia- 
tions from observations and synthetic spectra quantitatively, we 
have analysed the CCFs of r] Serpentis. The results are shown in 
Table|2]and shown in Fig.|9] For these stars we see again that m 
- -2 and m = 2 are the least favoured mode identification for all 
frequencies. 

For T] Serpentis all amplitude distributions are asymmetric 
and therefore none of the fits are consistent in the line wings, 
while in the centre both m - and m = -1 or m = 1 are con- 
sistent with the observed distribution. In the centre of the phase 
distribution the m - l(V and 2"'' frequency) and m = -1 (3"' 
frequency) fits match the observations closely, while the m - 
fits are just within the errors (1" and 2'"' frequency) or even out- 
side the error bars (3''' frequency). The phase distributions in the 
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Fig. 10. Same as Fig.|4]for 77 Serpentis. 

CCF wings are best matched by the m = fits. For this star the 
line centres favour an m = 1 ( 1 and 2'"' frequency) or m = - 1 
(3'''' frequency) mode identification, while m - cannot be ruled 
out because of the better fits in the CCF wings. 

The inclination angles of the synthetic line profiles fitted to 
the data have a weighted mean value of 57° with a standard devi- 
ation of 16° (see Fig.[TO]for the histogram). This range of incli- 
nation angles together with the projected rotational velocity de- 
termined with FAMIAS (Table |2]i results in a surface rotational 
frequency between 1.5 and 2.1 yuHz. 

6. Discussion and conclusions 

To take the analysis of line-profile variations performed by 
[Hekker et al. (2006) one step further, i.e., from a qualitative anal- 
ysis to a quantitative analysis, it was necessary to apply a number 
of specific corrections to the cross-correlation functions. With 
these corrections we removed profiles with outliers at any veloc- 
ity across the line profile. These outliers seem to have caused the 
large, single peaked amplitude profiles and low phase changes 
of the line-profile variations, most clearly seen in the dominant 
frequencies o f e Ophiuchi and 77 Serpentis, see Figs. 6 and 7 of 
lHekkeretan (2006). Hekker et al. (2006) could match these am- 
plitude profiles qualitatively with m - 2 modes, but the low phase 
change over the profile could not be reproduced from synthetic 
line profiles, which hampered a definite mode identification. The 
corrections applied here altered the amplitude and phase distri- 
butions and on these corrected distributions we could perform a 
quantitative analysis. 

We have been able to analyse line-profile variations in the 
corrected CCFs of four evolved stars and identify the azimuthal 
order for the frequencies detected in the first moments with an 
uncertainty of +1 as is usually the case in spectroscopic mode 
identification with the adapted method. From the synthetic line 
profile fitting of non-radial modes as performed with FAMIAS, 
we could also determine a range of inclination angles for the 
stars. These together with the projected rotational velocity pro- 
vided us with the surface rotational frequencies, which are all 
shown in Table [3] 

The uncertainty in the mode identification is mostly due to 
different favourable identifications in the line wings and in the 
line centres. Data with higher signal-to-noise ratio, i.e., lower 
noise levels, than we have at hand could improve the mode iden- 
tification significantly. Firstly, this would decrease the asymme- 
try in the observed amplitude distributions for modes with m - 



0, as seen from simulations. Secondly, lower noise values could 
reduce the errors on the amplitude and phase distributions which 
will allow us to better distinguish between fits with different m 
values. Nevertheless it is possible that higher signal-to-noise ra- 
tio of the data will not improve the fits to the phase distribution. 
It is clear from the lower panels of Figs. |5] Q and |9] that for m 
+ the phase difference in the wings changes gradually for the 
synthetic profiles, while this is less so for the observed profiles. 
We tested whether this discrepancy was due to the fact that we 
neglect temperature effects on the equivalent width, but taking 
those into account did not improve the fitting considerably. 

For j6 Hydri the azimuthal orders obtained from the line- 
profile analysis are compatible with the harmonic degree of 
the modes determined previously from radial velocity mea- 
surem ents and asymptotic frequency relations bv lBedding et al.l 
(I2007h . Three modes with f = and one mode with ^ = I are 
compatible with m = 0, while for the fifth frequency we know t 
-\ and ;« = -1 is clearly favourable in this case. 

The mode identificatio ns for 6 Eri d ani are also consis- 
tent with the results from Carried (l2002h : ICarrier et aP (l2003h . 
Although the present method seems to point to a non-radial 
mo de for the dominant frequency inst ead of a radial as claimed 
by Carried (l2002h : ICarrier et al.l(l2003h . Our identification would 
change the degrees assigned to the different ridges in the echelle 
diagram, but we do not have enough frequencies for which we 
could apply the current method to make any firm statements. 

For 6 Ophiuchi there is one mode in common between the 
present work and th e frequency model fitting performed by 
rKaninger 'etaD(l2008h . The identifications for this non-radial os- 
cillation mode are also in agreement. As we find two possible 
values of the azimuthal order for all frequencies, our analysis 
could also be in agreement w ith the radial mode interpretation 
favoured bv'De Ri dder et all (12006) and'B arban eTall ([2007). 

The confirmation that we can compare line-profile variations 
of observations and synthetic profiles quantitatively led us to 
also analyse r\ Serpentis. Also for this stars we were able to ob- 
tain the azimuthal orders for three modes with an uncertainty of 
+ 1 and an indication of the surface rotational frequency. These 
identifications and the surface rotational frequency are less accu- 
rate due to lower signal to noise in the amplitude and phase di- 
agrams which are therefore more asymmetric than for the other 
stars. 

With FAMIAS we fitted usin/ independent from known 
literature values. For all four stars we find higher values than 
current literature values, see Tables [1] and |2] The difference 
in the quoted values might be due to differences in equivalent 
width and^niacio- 

As discussed here, improvements in the spectral line profile 
analysis are still needed, both in terms of higher signal-to-noise 
ratio observations and in terms of the generated synthetic 
profiles we compare the data with. Expandi ng on an analy sis of 
simulated data, as akeady started by Hek ker et al.l (l2006h . and 
applying the method to evolved stars with solar-like oscillations 
observed with high signal-to-noise ratio as we expect to be- 
come available from for instance SONG (Stellar Observations 
Network Group), will improve our understanding and increase 
the value of this method for mode identification in solar-like 
oscillators. 
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Fig. 9. Same as Fig.|5]for the significant frequencies of rj Serpentis obtained from (v): V(v> = 11.74 cd ' (135.9 //Hz) (left), V(v> = 
9.48 cd"' (109.7 yuHz) (centre) and v^v) = 10.38 cd"' (120.1 yuHz) (right). The mean radial velocity of the star is approximately 
9.4 km s"'. 



Table 3. Intervals for the surface rotation frequencies (Q) for 
the four evolved stars discussed here, together with the incli- 
nation angle (/) and projected rotational velocity (v sin /), we 
used to compute Q.. 



star 




n 






u sin ; 






de 


jree 


kms"' 


j3 Hydri 


3.6 


-5.5 


38 


-72 


4.3 


e Ophiuchi 


0.8 


- 1.1 


45 


-81 


5.7 


Tj Serpentis 


1.5 


-2.1 


41 


-73 


5.6 


S Eridani 


3.1 


-4.9 


38 


-76 


4.9 
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